Separation of spacetime and matter in polar oscillations of compact stars 
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In polar oscillations of compact stars dynamical changes in spacetime and matter are in general 
coupled together. Based on the theory of stellar pulsation proposed by Allen et al. [Phys. Rev. D 
58, 124012 (1998)], we establish a feasible approximation scheme to separate spacetime from matter 
in such oscillations. Two independent second-order ordinary differential equations are obtained, 
which lead to accurate determination of p-mode and UMnode quasi-normal modes of compact stars, 
respectively. 
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CN \ I. INTRODUCTION 

bTj 

Compact stars, including neutron stars (NSs) and quark stars (QSs), are notable for the high density achievable 
at their centers and the huge gravitational field surrounding them. The static structure of NSs (or QSs) is often the 
ideal test bed for the equation of state (EOS) of nuclear (or quark) matter and can provide clues to the interactions 
between elementary particles as well (see, e.g., [H-Q and references therein). Likewise, gravitational waves (GWs) 
emitted from oscillating NSs (or QSs) are also expected to carry rich information about the internal structure of these 
stars, giving rise to various efforts to infer the properties of relativistic compact stars from close examination of their 
oscillation spectra (see, e.g., 0412 )• Although the oscillation frequencies of NSs usually lie outside the sensitivity 
window of the existing GW detectors such as GEO600, LIGO I, TAMA300 and VIRGO (see, e.g., (HQ), with thc 
£^ advent of more sensitive GW telescopes (e.g., the proposed Einstein Telescope [H, Hp), it might be possible to infer 
i— —i. the structure and the composition of a compact star from its GW signals in the coming future. 

The fully relativistic formulation of stellar pulsations was founded by Thorne and Campolattaro in 1967 [ij} and 
since then there have been other attempts to describe such processes with other alternative approaches (see, e.g., 
[l8l - |22j ] ) . Basically, oscillations of compact stars are classified into axial and polar types. For axial type oscillations, 
\ the matter forming a star is essentially a spectator and only the spacetime surrounding the star pulsates in accordance 
with a scalar wave equation [2Cj| . On the other hand, in polar type oscillations matter motion and spacetime variation 
are coupled together and hence polar type oscillations are expected to be more relevant to dynamical processes such 
. as supernova and binary mergers of NSs. 

To describe the interplay between matter and spacetime in polar oscillations of compact stars, different versions 
£NJ , of equations of motion governing such cou pled motion have been proposed since the pioneering work of Thorne and 
t— \ • Campolattaro [17j. For example, based on [17[, Lindblom and Detweiler described the system with four coupled first- 
order differential equations (termed as the LD formalism in this paper) [H|,[l9j], which can be readily solved numerically. 
• i-h ■ Later, to draw an analogy between the oscillations of relativistic and Newtonian stars, Lindblom, Mendell and Ipser 
(2l| proposed a new formalism to describe the coupled oscillation with two second-order equations (termed as the 
LMI formalism in this paper) , which reduce directly into the corresponding Newtonian equations under the weak field 
- - " limit. 

In another relativistic stellar pulsation formalism proposed by Allen, Anderson, Kokkotas and Schutz (termed as 
the AAKS formalism in this paper) [23], variations in matter and spacetime are described by one matter variable 
and two metric variables, which are governed by three coupled second-order time-dependent wave equations. In fact, 
these three variables are also related through a time-independent Hamiltonian constraint. Hence, only two of the 
three wave equations are independent. Allen et al. [23| managed to evolve these variables in the time domain. 

The aim of the present paper is to study polar pulsations of compact stars with the AAKS formalism in the 
frequency domain where all dynamical variables are assumed to have exp(— iuit) time dependence, with lo being the 
eigenfrequency. Oscillation modes in which GWs are emitted from a star to spatial infinity are termed as quasi-normal 
modes (QNMs) because they are usually damped [2~4T - |28l ] . Hence, the eigenfrequency u) = u> T — iuii is complex- valued 
with Wi > measuring the decay rate. QNMs of compact stars are worth studying because, as mentioned above, they 
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can reveal the physical characteristics (e.g., mass, radius, moment of inertia, equation of state (EOS) and composition) 
of a star (see, e.g., (3-TTT| ) . On the other hand, polar QNMs of compact stars can be classified into the fluid mode 
(including the fundamental /-mode, the pressure p-mode and the gravity (/-mode according to the nature of respective 
restoring forces (29j |) and the spacetime iw-mode which has no Newtonian counterpart. 

In the present paper, based on the AAKS formalism, we derive in the frequency domain two coupled wave equations 
governing the evolution of the matter field and the metric variables. As in the LD formalism, these two equations 
together lead to both fluid and spacetime QNMs of compact stars. In addition, we show that these two equations 
can be straightforwardly decoupled, resulting in two independent equations which can reproduce good approximation 
of fluid modes and w-modes, respectively. With such decoupling scheme for matter and spacetime, we succeed 
in obtaining the Cowling approximation (CA), where effects of metric perturbations are neglected HH, and the 
inverse-Cowling approximation (ICA), where fluid motion is ignored [3ll.l32j. for compact stars in a unified framework. 
The fluid motion is described by a single second-order equation in the matter field without involving any spacetime 
perturbations. Meanwhile, polar w-mode oscillations obey a single wave equation as in the case of axial w-modes [2~oj]. 
To our knowledge, this is the first time a single second-order equation for polar w-mode QNMs is found. 

Furthermore, we show that the decoupled equation for fluid motion, upon imposing suitable physical boundary 
conditions, form a Sturm-Liouville eigenvalue problem (see, e.g., [33jj ) . Therefore, the fluid system can be studied with 
conventional normal-mode analysis, and the energy of the fluid is conserved as long as the influence of the spacetime 
perturbations are negligible. The mode with the lowest eigenfrequency corresponds to the /-mode oscillation, while 
other modes are good approximation for p-mode QNMs. 

The organization of the present paper is as follows. We briefly review the equilibrium state and the perturbed 
state of a compact star in Sees. II and III, respectively. The AAKS formalism for pulsation of compact stars in the 
frequency domain is developed in Sec. IV. In Sees. V and VI we recast the equations of motion in the AAKS formalism 
into two coupled wave equations in two variables. We show in Sees. VII and VIII that these two equations can be 
decoupled, respectively leading to accurate CA and ICA schemes. We conclude the present paper with a discussion 
in Sec. IX. Besides, there are two appendices respectively providing detailed discussion on the boundary condition 
imposed at the stellar surface and an alternative approach to description of pulsation of compact stars. 

Geometric units in which G = c = 1 are adopted throughout this paper. We consider polar oscillations of a non- 
rotating compact star with radius R and total mass M. Unless stated otherwise, all numerical values of frequencies 
are evaluated for the case where the angular momentum index I equals 2 and measured in units of M . Besides, we 
also assume that the effect of temperature on the EOS of nuclear matter is negligible and hence g-modes are omitted 
in our discussion. 



II. EQUILIBRIUM STATE OF COMPACT STARS 



The energy-momentum tensor of a non-rotating compact star made of a perfect fluid at equilibrium state is given 

by 

T Q/3 = {p + p)u a u p +pg a ^ (2.1) 

where p is the pressure, p is the energy density and u a is the 4- velocity of the fluid. The spacetime around the star 
is described by the line element: 

ds 2 = g a(s dx a dx = -e u dt 2 + e x dr 2 + r 2 (d9 2 + sin 2 ed(f> 2 ), (2.2) 

where (x ) = (t,r,6,(/>), and reduces to the standard Schwarzschild form outside the star, 

ds 2 = -(1 - —)dt 2 + (1 - —y l dr 2 + r 2 (d9 2 + sin 2 9dq> 2 ). (2.3) 
r r 

By substituting T Q/3 and g a p into the Einstein equation, the Tolman-Oppenheimer-Volkov (TOV) equations of 
hydrostatic equilibrium can be derived [H], IH| : 

dp _ (p + p)(m + 4?rr 3 p) ^ 



dr r(r — 2m) 

dv 2(m + 47rr 3 p) 



(2.5) 



dr r(r — 2m) 

e~ x = 1 - 2m/r, (2.6) 
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with 



m(r) 



47rr' 2 p(r')dr' 



(2.7) 



The hydrostatic structure of a compact star can be specified by the TOV equations, together with the star's EOS 
P = p(p) an d the boundary conditions: 



e v \ r =R = 1-2M/R, 

p\r=R = 0. 



(2.8) 
(2.9) 



III. PERTURBED NEUTRON STARS 

When a static compact star is perturbed, its fluid and spacetime perform small-amplitude oscillations about the 
equilibrium state. Hence, the perturbed metric is given by: 



Ja/3 



(3.1) 



where g a p is the spacetime metric of the equilibrium configuration given by the TOV equations [3J, |35[ . The per- 
turbations in the spacetime metric and fluid displacement are first decomposed into tensorial and vectorial spherical 
harmonics with angular momenta {lm}, and then separated into axial and polar parts accordin g to parity. The polar 
part of the perturbation in the metric, hf , in the Regge- Wheeler gauge can be written as [l7l.l36l.l37j: 



(him) al- 



and the corresponding fluid displacements are |17Ml9j : 



/ e u H H x 
Hi e x H a 0_ 
r 2 K 



Or 2 sin 2 OK J 



e = ~r- 2 Vd e Y lm , 

^ = -{rsme)~ 2 Vd^Y lm , 



where Hq, Hi, K, W and V are functions of both t and r. 



IV. AAKS FORMALISM IN FREQUENCY DOMAIN 
A. Equations inside the Star 



(3.2) 



(3.3) 
(3.4) 
(3.5) 



In the AAKS formalism [23j, polar oscillations of compact stars are described by two spacetime metric variables 
F(t, r), S(t, r) and one fluid perturbation variable H (t, r) which are related to the metric and the Eulerian change in 
pressure, dp, as follows: 



S(r,t) 



e v [H (r,t) - K{r,t)} 



F(r,t) = rK(r,t), 
Sp(r, t) 



H(r,t) = 



P + P 



For QNMs, the three variables S(r,t), F{r,t) and H{r,t) all have the e lu)t time-dependence, i.e., 

S(r,t) = S(r)e-^, 
F(r,t) = F{r)e-' lut , 
H{r,t) = H{r)e- Wt . 



(4.1) 
(4.2) 
(4.3) 



(4.4) 
(4.5) 
(4.6) 
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It is readily shown that S(r), F(r) and H(r) satisfy the following three coupled second-order differential equations 
231: 



d 2 S 2e v 

S+~^j + — [2nr 3 {p + 3p) + m- (n + l)r] 5 



ie 



2v 



3\2 



(m + Anpr ) 
r — 2m 



47rpr — 3m 



i 7 . 



d 2 .F 2e l/ 



> 2 F + — + — [27rr 3 (3p + p) + m - (n + l)r] F 
-2 [47rr 2 (p + p) - e~ A ] S + 8n(p + p)re" (1 - -^)H, 



d 2 H e^+ A )/ 2 dii 



dr? 



r 2 dr* 



27rr 2 (p + p)(3 + ^)-(n+l 



if 



(m + Anpr 3 ) 2 , 1 . m + Airpr 3 , 
— — (1 H ) — (1 • 

r 2 (r-2m) 1 C 2 ' 2r 2 K 

2(m + iirpr 3 ) 2 1 m + 4-irpr 3 

^ (1 



r 2 (r-2m) C 2 



(m + 47rpr 3 )(l - i) + 2(r - 2m) 

) — 47rr(3p + p) 
) — 47rr(3p + p) 



2r 2 



1 

Cf 
1 

C 2 



s 

F 



+(m + 4Trpr 3 )(l- — ) 



1 e ( A -")/ 2 / e " dF d5\ 



C 2 2r V r 2 dr* dr* 7 



where 



n= (Z + 2)(I-l)/2, 

r* is the standard tortoise coordinate related to the circumferential radius r through the relationship: 



r*{r) 





'X(r')-u(r'Y 


/ exp 
Jo 


2 



dr' 



and C 2 (r) is the square of the sound speed in the stellar fluid, defined as 

C 2 = dp/ dp. 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.H) 



(4.12) 



In addition to the above three equations, S(r), F(r) and H(r) should also satisfy the following Hamiltonian 
constraint: 

d 2 F e (l/+A)/2 dF e v 

— 7 (m + 4irr 3 p)- h —[12irr 3 p - m - 2{n + l)r]F 



dr 2 
rp C+ x )/ 2 



dS 
dr* 



inr 2 (p + p) - (n + 3) + — 

r 



S + ^e»(p + p)H = 0. 



(4.13) 



By combining (|4.8[) and the Hamiltonian constraint (|4.13[) together, a simple expression for H can be derived |23j |: 

^ {X+U)/2 - dF „ . e",_ . . ,„^,o dS 



87t(/j + p)re u H 



(m + Airpr 3 )— + w 2 F + — (3m + ^pr 3 )F + e -^+ x )/ 2 r — + (n + 1)5, (4.14) 
dr 



which is often used to calculate the value of H, instead of integrating (14.91) . and termed as the modified Hamiltonian 
constraint hereafter. However, Eq. (|4.14|) obviously fails to give the value of H at the stellar surface, where both p 
and p vanish. We will further discuss such subtlety in the later part of this section. 
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B. Equations outside the Star 

As for the region outside the star the fluid variable H vanishes identically, the two metric variables S and F evolve 
as: 

^ S + + 2(r ~ 4 2M) [M - {n + l)r] S = -^(r - 2M)(7M - 3r)F, 

(4.15) 

W 2 F + + [M - („ + 1),] F = 2(1— ™)S, (4.16) 

which can be obtained from (|4.T[) and (|4.8|) by setting p = p = and m(r) = M. In addition, they arc related to the 
Zerilli function Z(r*) as follows: 

(n + l)(»r + 3M) z = _ 1 dF + (n + 2)r M p + 
r 6 r dr* r A 

where the Zerilli function Z(r*) satisfies the Zerilli equation (38j : 

d 2 



Z(r*) = 0, (4.18) 



with the Zerilli potential Vz(r*) given by 



r — 2M 

Vz(n) = -J7 \2n 2 (n + l)r 3 + 6n 2 Mr 2 + !8nM 2 r + 18M 3 } . (4.19) 

r 4 (nr + &My 

C. Evaluation of quasi-normal modes 

QNMs of vibrating compact stars arc defined by the physical solutions of the AAKS equations with the Zerilli 
function satisfying the outgoing boundary condition Z — > exp(iwr*) at spatial infinity. In solving the AAKS equation 
set inside the star, we numerically integrate (|4.7|) . (|4.8|) outward to find S and F, and use the modified Hamiltonian 
constraint (|4.14p to evaluate H. However, a few remarks are noteworthy in carrying out the numerical integration 
scheme from r = to r = R. 

First of all, there is a regular singular point at r = in (|4. T[) - (|4.9[) resulting from the term (n + l)/r 2 , which 
forbids normal discretization scheme. Instead, the regularity boundary condition for S, F and H at r = leads to 
the following asymptotic behavior near the origin: 

S « s r l +\ (4.20) 
F a f r l+ \ (4.21) 
H a h Q r l , (4.22) 

where f , s and h are constants. Expanding the modified Hamiltonian constraint (14.141) about the origin, we find 
that these three constants satisfy the following relationship: 

2(po+po)e uo h = [e Uo/2 (po/3+p )(l + l)+uJ 2 +e uo (po+po)}fo + [e- uo/2 (l + l)+n + l}s , (4.23) 

where po, po and e~ v ° stand for the values of p, p and e~ v evaluated at r — 0, respectively. We make use of the 
asymptotic behavior in (|4.20[) . (|4.21j) and (|4.22j) and start numerical integration at a point slightly away from the 
point r = 0. Owing to the relationship (|4.23p . two independent solutions can be obtained when the AAKS equations 
are integrated outward. We then impose a proper boundary condition at the stellar surface r = R to fix the ratio of 
these two independent solutions. 

Since the Lagrangian perturbation of pressure Ap vanishes at the stellar surface, we have 

Sp(r = R) = - C^r ■ ( 4 - 24 ) 

L dr \r=R 
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The value of H there then follows directly from the definition of H, (|4.3I) . and the TOV equations [391 ]: 

Together with the modified Hamiltonian constraint (14. 14|) and the fact that p = p = at the stellar surface, the 
boundedness of H (r = R) implied by the above equation leads to the following boundary condition at r = R: 

However, we find that this is merely a necessary condition, instead of a sufficient one, and cannot be considered as 
an independent boundary condition. 

Instead, another independent boundary condition can be obtained by keeping only terms of order in (|4.9|) as 

C s — > when approaching the stellar surface [39|, |4CJ : 

MdH M dF M dS M(R + 2M) M 

R 2 dn 2i? 3 dr* 2{R - 2M) dr* 2i? 5 2R(R - 2M) { ' 

By imposing this condition on the linear combination of the two independent solutions, we obtain a single (up to a 
multiplicative constant) solution of S, F and H inside the star, which is connected to the Zerilli function outside the 
star to determine the QNM frequencies. 

However, a remark about the accuracy of the direct integration scheme mentioned above is in order. Since the 
modified Hamiltonian constraint (|4.14l) is used to find H, considerable numerical errors could arise near the stellar 
surface where both p and p are small. To overcome this problem, we analytically expand S, F and H as power series 
there and connect such expansion with the solution obtained from numerical integration starting from the center of 
the star at a point near the stellar surface (see Appendix [A} . We find that such a matching scheme can yield more 
accurate results than the direct integration scheme. Moreover, the expansion could shed light on the relationships 
among the three physical fields (S, F and H) near the stellar surface, which are useful to the formulation of the 
decoupling schemes developed in the later part of the present paper. 



V. H-S APPROACH 



As mentioned in the Introduction, polar QNMs can be divided into two categories, namely, fluid modes and 
spacetime modes. In fluid modes, especially the p-modes, the associated GW emission is almost negligible, while in 
the spacetime w-modes, the fluid is hardly involved. It strongly suggests that appropriate schemes can be developed 
to decouple spacetime oscillations and fluid motions, leading to reduction in the order and hence simplification of the 
system. 

It has been shown that inside the star and in the weak field limit K and Hq are related to the perturbation in the 
Newtonian potential, U\, by |4lj |: 

K ~ iJ ~ -2U X . (5.1) 

Therefore, following directly from (|4.1[) and (|4.2[) . to leading orders in field strengths, S and F measure the relativistic 
and Newtonian effects of gravity, respectively. Hence, in order to describe generation of gravitational radiation, which 
is the dominant process in the spacetime mode, S is the right variable to keep. On the other hand, H is proportional 
to the Eulerian change in pressure and describes the motion of matter. 

In the following, based on the AAKS formalism, we establish a two-equation approach and a decoupling scheme 
leading to accurate CA for the fluid modes (especially the p-modes) and ICA for the iyj-modes using H and S, which 
is termed as the H-S approach in the present paper. First of all, we express F in terms of S and H for the regions 
r < R and r > R. 

For the case r < R we rewrite (|4. 14|) as: 

6n(j> + p)e^rH = A f ~ + A f0 F + A sl ^- + A s0 S, (5.2) 



with 



A sl = e-^/\ 
Ao = n+ 1, 

e (A+^)/2 

A/i = — — (m + Anpr 3 ), 



Af Q = uj 2 + -3 (3m + Airpr 3 ). 



Differentiating (|5.2I) once, and using (|4.7[) and (|4.8|) to eliminate d 2 F/dr 2 and d 2 S/dr 2 , we get: 



dF 



rfS 



8tt( P + p )e» — {rH) = B fl — + B fQ F + B sl — + B s0 S, 



where 



B s i = 1 + n + e A - 4?rr 2 (p + p), 

e (^-A)/2 

B s0 = [2(1 + n) -e^V^ 2 - 47rr 2 (p + p) (1 + rV)] , 



/1 



2 1 1/ 

cj + e 



(«/-A)/2 



/0 



4?r(2p + p) 



e w (3m - 47rr 3 /9) ^~ ~ ^"t 1 + n ~ 87rr ' 2 (P + P))] 



with v' — dv J dr. 

From (|5.2p and (|5.7I) . we can express F and dF/dr* in terms of S 1 , dS/dr*, H and dH/dr*: 

dS a dH 

r = a sl - h a s0 b + a h i- \-a>hoB, 

dr* dr* 

dF a dS c a dH a rr 

-7- = Psi-r- + PsoS + 0hi^— +PmH, 
dr* dr* dr* 



where 



or.l = -^\er^^ 2 ru 2 + \e (u - x ^ 2 \4nr 3 ( P + p)-e x (m + 4TTr 3 P )(l + n) 
D r z 



+4iTr 2 e x (to + 4nr 3 p) (p + p)] >, 



as0 = ~D 



n + 1 H h 47rr p 



+ — e v [-(1 + n)(m - 4nr 3 p) 



Othl 



+47rr 2 e A (m + 47rr 3 p)(l + 87rr 2 p)(p + p)] ^, 
8 ^ e (A+3,)/2 (/g + p) 



-(to + 47175?" ), 



87re !/ (p + p)r f 2 -, 
- |w + 47re [p + p)\ 



D 



Pel = — R 



n + 2 4-7rr 2 p w 2 + — e" {m [l + n - 4irr 2 (p + p)] 

+47rr 3 [p-np+ \2i:r 2 p(p + p)] } j, 

_ _L e -(A+^)/2 J _ rw 4 + 1 e „ | g A/ m + 47rr ,3 p ) N + n _ 87r? ,2/ + 

D y r 

1 f 2m 

+2r(n + 1) - 3m - 47rr 3 (2p + p)} w 2 + —e 2v I mfl + n H ) 

+2(m + 47rr 3 p) — + n-e A (l + n) 2 + 47rr(p + p)(m + 47rrV) 
L r J 

{r - 2e x [m - r(l + 2n) + 47rr 3 p] } + 47rr 3 p[l 



(5.18) 



2m. 



w 2 + — (3m + 47rpr 3 ) 



87re (3^A)/2( p + p ) 



1 + -e A (m + 47rr 3 p) 



• 2 + e l/ -^{4 7 rr 3 (p + p) 



-2e x (m + 4?rr 3 p) [l +n - 87rr 2 (p + p)] } L 



(5.19) 
(5.20) 

(5.21) 



and 



2e 2 



■ {87rr 3 (m + 2irr 3 p)(p + p) - e A (m + 47rr 3 p) 2 [l + n - 87rr 2 (p + p)] } 



+— [e A (m + 4nr 3 p) 2 + Amr + 4?rr 4 (3p + p)] w 2 



- w 



(5.22) 



For r > R, the stellar pressure p and density p vanish and so does H, then F 1 in (|5.12j) reduces to the following 
form: 



(c) d5 (e ) „ 



(5.23) 



which is a linear combination of 5 and dS/dr*. The coefficients and are obtained from a s o and a s i by taking 
p = p = and m = M: 



(e) 

aw = — 



1 



1 



(e) _ 

tt s0 - D(e) 



(n+l + -Af)w 2 - — M(l + n)(r - 2M) 



with I)( e ) being: 



D [c) = oj 4 + —M(-7M + Ar)u 2 - 4rA/ 2 (l + n)(r - 2M). 



(5.24) 
(5.25) 

(5.26) 



VI. EQUATIONS OF MOTION IN H-S APPROACH 



With the expressions of F and dF/dr* derived above, we manage to rewrite the AAKS formalism in terms of the 
fields H and S, which is useful in decoupling matter and spacetime in oscillations of compact stars. 



A. Equations inside the Star 



The system inside the star can be described by two coupled second-order ODEs, 



d 2 S rS dS 



drl 
d 2 H 
dr 2 



? ~ Vf - Vf a s0 )S = V/ (a u ^ + a h0 H), 



dH 



dH 

dr* 



+ (Vftl - V? x fax - V>hx)^r + (^2 - V& - Vftpho - V f H a h0 )H 



dS 



= (Kf + Vfxfa + y^ sl )— + (V% + Vflfa + V f H a s0 )S, 
where the a and f3 coefficients are defined in (|5.14[) - (|5.21[) . and 

V% = %■ [(» + l)r - 2nr 3 (p + 3p) - m] , 



V s 



v hl — 

v H - 

v h0 — 



4e 



2;/ 



(to + Airpr ) 



3\2 



,(f+A)/2 



2to 



+ 47rpr — 3to 



1 



2e v 



(to + 47rpr 3 )(l - — j ) + 2(r - 2m) 



l-2^r 2 ( p + p)(3+— ) 



1 e ( A+,y )/ 2 
y^ = ( m + 47rpr 3 )(l-— ) 



V f0 



C s 2 y 2r 3 ' 
2(TO + 47rpr 3 ) 2 1 to + Airpr 3 



r 2 (r — 2m) C 2 



2r 2 



(to + 47rj?r 3 )(l 



1 e (A-v)/2 



c 2 



(m + iirpr 3 ) 2 , 1 , 
- — —(1 H ) 

r 2 (r~2m) V C 2> 



2r ' 
to + Anpr 3 
2^2 



(1 - ^2 ) - 47rr(3p + p). 



The S'-equation (|6.ip can be recast into a much nicer form with the following transformation: 

S(oj, r») = <?(w, r*)5(w J r*), 

where q(uj,r*) is defined as: 



= exp 



Vf (ri)a,i( W y»)dri 



In terms of the transformed variable S(uj,r*), Eq. ()6. 1[) becomes: 



d 2 S 



— + [ w 2 -y« (Wjn)] 5 = ^ 



with 



V«(u;,r,) = Ko(^*) + Vf (r.)a M (u,r.) 



1 d 

2 dr* 



[Vf (r*)a s x(u,r*)] + - [V^(r*)a s i(o;,r*)] , 



Oh 



:Vfo(n) 



a w (o;,r*)— + a/jo(u;,r*) 



s ?(^,r*) 

The ff-equation (|6.2p can also be written in terms of S, 



^ + (V h H x + 5V h »)f- + - ^ - 5V&)H = H S, 



(6.1) 
(6.2) 

(6.3) 
(6.4) 

(6.5) 

(6.6) 

(6.7) 

(6.8) 

(6.9) 

(6.10) 

(6.11) 
(6.12) 

(6.13) 

(6.14) 
(6.15) 

(6.16) 
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with 



SVhi = -V f H tl3hi-V f H a hl , 
= V f H ^ h0 + V f H a h0 , 

Oh = (Kf + Vj[0 al +V%a. 1 ) 



dq(uj,r*) + ^ d 
dr* ' * dr* 



(6.17) 
(6.18) 

(Ko + VfiPsO + V f H a s0 ) q(u, r*). (6.19) 



B. Equations outside the Star 

Since the functions a^i(w,r*) and aho{uJ,r*) both contain a (p + p) factor, they vanish outside the star. Therefore, 
the system reduces into one second-order ODE in S outside the star: 



rS(c) T/ S(e) (e)x Q _ . 



where 



K 5 o (o) = ^(r-2M)[(l + n)r-M], 
V^ e ) = _ (r-2M)M(3r-7M), 



and the c^ ' coefficients are defined in (|5.24[) and (|5.25|) . 
As a counterpart of (|6.13[) . S satisfies 



(6.20) 

(6.21) 
(6.22) 



d 2 S 
dr 2 



+ [lo 2 -V( Q) (w,r,)]S = 0, 



where 



3.23) 



= 2( ^~^y |2M 4 (n + l) 2 (2Af - r)[3M 2 - 2r 2 (n + 1) + 2Mr(2n + 1)] 

-M 3 r 3 (n + 1)cj 2 [126M 3 - 246M 2 r - Mr 2 (8n - 131) + 4r 3 (n - 5)] 
-Af 2 r 6 w 4 [147Af 3 - M 2 r(97n + 258) - 4Mr 2 (2n 2 - 17n - 33) + r 3 (4n 2 - lln - 21)] 

-Mr w uj 6 [56M 2 - hlMr - 2r 2 {n - 5)] + r 14 w 8 [(n + l)r - M] 



3.24) 



and 5 is related to S through 



with g( e )(w, r*) given by: 



q (c) (w, r*) = <?(w, exp 



3.25) 



(6.26) 



Hereafter i?* is the tortoise coordinate of the stellar surface r — R. From the forms of Vj?^ and , it is easy to 
see that approaches a constant at the spatial infinity. 

It is interesting to note that the Zerilli function is expressible solely in terms of S (or equivalently S) in the H-S 
approach as H is identically zero outside the star, which reads as: 



( f 
\ dZ/dr* 



Zii(w, r) Zu(w, r) \ / S 
Z 2 i(uj,r) Z 2 2(u,r) J \ dS/dr* 



3.27) 
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with 



Zu 

Z2I 



-22 



1 f M(r - 2M) [9M + (-3 + n) r] 



Z)( c ) \ (3M + nr)r 4 

1 M[3M - (3 + n)r] 
Z)( c ) (3M + nr)r 2 ' 
M [-9M 2 + 9Mr + n(n + 3)r 2 ] u? 
D(c)r 2 (3M + nr) 2 
M(r - 2M) 



■ or 



D( e )r 6 (3M + w) 2 



[54M 3 + 3(14n - 3)M 2 r + 6n(3n + l)Afr 2 + n(n + l)(2n - 3)r 3 ] , 



3M(r - 2M) [3M 2 + InMr + n(n + l)r 2 ] 
r 4 (3M + nr) 2 



(6.28) 
(6.29) 



(6.30) 
(6.31) 



Therefore, Eq. (|6.20l) or Eq. (|6 . 23[) can be considered as an appropriate replacement for the Zerilli equation in the 
current theory. Likewise, S or S can in fact play the role of the Zerilli function Z. For example, QNMs can be defined 
by the imposing the outgoing wave boundary condition on S (or S). 

To sum up, in the H-S approach of the AAKS formalism a pulsating NS is governed by two coupled second-order 
equations, (|6 . 1 3|) and (16.16[) . inside the star, and one single second-order equation (|6 .23[) outside. The polar QNMs 
can be located by imposing the regularity boundary condition at the star center and the outgoing boundary condition 
at the spatial infinity. In the following discussion, we will show that the H-S approach achieved here forms the basis 
of a decoupling scheme readily leading to accurate CA and ICA for p-modes and w-modes, respectively. 



VII. COWLING APPROXIMATION 



The p-mode pulsation is one important type of fluid modes, whose restoring force comes mainly from pressure, and 
has been well studied for non-relativistic stars under the CA (29|. Here, based on the H-S approach developed above, 
we propose a relativistic version of CA for p-modes of compact stars. We will also see that such a CA scheme could 
also provide qualitatively correct description for the /-mode 



A. Fluid equation 



It is well known that the eigenfrcqucncics of the p-modes of a compact star are characterized by dramatically small 
imaginary parts. The smallness of the imaginary part implies the tininess of GWs emitted in p-modes. Therefore, 
the influence of the variable S on the fluid motion is expected to be negligible when calculating p-mode frequencies. 
With this in mind, we propose that the approximate eigenfrequencies of p-modes can be calculated with the following 
equation: 

+ (Vhi + SV&)%- + - V& - SV h " )H = 0, (7.1) 

which is obtained by neglecting the S part in (|6.16|) . Furthermore, we find that for p-mode oscillations the effects of 
SVu[ and SV^ are pretty small. Then, Eq. (|7.1I) can be further simplified into: 

The above equation can also be derived by simply neglecting the S, dS/dr*, F and dF/dr* terms in (I4.9[) . Unlike 8V& 
and 8VfH, which are frequency-dependent, V^q and V^f are both frequency-independent. Notwithstanding the fact 
that Eq. (|7.1I) properly accounts for the effects of F on the fluid motion and is proved to be the best approximation 
for p-modes in the following discussion, the mathematical structure of (|T. 2[) is more desirable and leads to a complete 
orthonormal set of eigenfunctions. 
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B. Boundary conditions 



Under CA (|7.2p . p-mode oscillations are determined from (|7.2p by imposing proper boundary conditions on H at 
r = and r = R. Following directly from (|4.22p . H w /lor", near the star center. Assuming that the EOS acquires the 
polytropic form p oc p 1+1 / N with the polytropic index N > 1 near the star surface, we show that (see Appendix [AT) 

C 2 S ~A C (R* -r„), (7.3) 

with ^4 C = M/(R 2 N). Consequently, the iJ-equation near the stellar surface reduces asymptotically to 

d 2 H M dH oj 2 

dr 2 R 2 A c {R*-r*)dr* + A C (R* - r*) ' ( ' j 

There are two independent solutions for this equation. The first one is bounded and behaves as: 

H oc [l + w 2 i? 2 (r* -R*)/M] . (7.5) 

The other one is proportional to (i?* — r*) 1_Ar and is unbounded (also see Appendix |A|) . As argued previously, H is 
finite at the stellar surface and therefore follows (|7.5I) . 

The single ii-equation (IT.2|) together with two boundary conditions (|4. 22[) and (|7.5|l form an eigenvalue problem. 
The approximate eigenfrequencies of polar p-modes under the CA can then be located. 



C. Completeness and orthogonality of p- modes 

The eigenvalue problem of the purely matter-based -H-equation (|7.2p is actually a Sturm-Liouville eigenvalue prob- 
lem, as it can be cast into the following standard form: 

-Q(r*)ff = -cj 2 A(r*)ff, (7.6) 

with 

P(r t ) = e v (p+p)r 2 , (7.7) 

Q(r.) = 2e 2 »(p + p)[n+l-2nr 2 (p + p)(3 + l/C 2 )], (7.8) 

A(n) = e»(p+p)r 2 /C 2 . (7.9) 

There are a series of real eigenfrequencies {<t„} and the corresponding normalized eigenfunctions {i?„(r !t )}, which are 
defined by the boundary conditions derived above, namely (i) H is regular at r* = 0, (ii) (dH/dr*)/H = R 2 ui 2 /M 
at r = i?*. Following straightforwardly from the standard theory of Sturm-Liouville system and the fact that 
A(r = 0) = = A(r = R), the normalized eigenfunctions {-£/«(?%)} form a complete orthogonal set obeying the 
completeness and orthogonality relationships (see, e.g., [33l|): 

tyv-O; (7-10) 

*nm- (7-11) 



P(r.) 



^jA(r*)ff„(r*)if m (rl) = 




dr*A(r*)H n (r*)H m (r*) = 



D. Numerical results and discussions 



To gauge the accuracy of the CA developed (or ICA as discussed later) in the present paper, we use EOS A [42[ 
to construct a NS with compactness C — 0.27 and central density p c — 2.227 x 10 -3 and evaluate the QNMs of 
the star. Table fl] shows the frequencies of the leading fluid modes obtained from exact calculation and CA (|7.2[) 
developed above. We see that the lowest eigenfrequency a of (|7.6j) provides us an approximate value of the real 
part of /-mode frequency, w r j, with a percentage error of around 20%. For n > 1, the eigenfrequency a n yields an 
accurate approximation for the real part of the frequency of the n-th p-mode, w r ,n- The percentage error is usually 
much less than 1% and decreases with increasing n (see Table U. 
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The comparison between 77-function obtained from (17.2[) and the exact ii-function is shown in Fig. [TJ Both of 
them are normalized so that 

R. 

dr»A(r„)|J3f = 1. (7.12) 

Besides, we assume that the exact i7-function is real at r = R. The numerical results again confirm that the CA 
proposed here is very accurate for p-modes. However, the discrepancy between the approximate and the exact wave 
functions grows larger for /-mode, which implies that the coupling between spacetime (i.e., F and S) and matter is 
negligible for p- modes, but is not so weak for /-mode. Actually, from the form of (|6.2I) . it is clearly seen that the 
accuracy of the CA (|7.2[) . would decrease if D(r,u>) is small. As shown in Fig. [5J where D is plotted against r/R 
for different values of ui, D decreases with decreasing lu. In fact, for u = <r (i.e., the /-mode frequency under CA), 
D vanishes at r/R = 0.947. Therefore, the omission of the influence of spacetime (i.e., F and S) on oscillations of 
matter field may lead to appreciable errors at low frequencies. On the other hand, we observe that there are exactly 
n nodes (excluding the origin) in the wave function H n , in agreement with standard theory of the Sturm-Liouville 
eigen-system (see, e.g., |43j|). 

We also note that the high accuracy of the CA (|7.2p is generic and independent of the EOS and the compactness 
of the NS. Tables [TTI and Hill show the exact complex eigenfrequencies (oj t ,lj{) and approximate real eigenfrequencies 
a„ obtained from CA scheme (|7.2p for the leading fluid modes of NSs with Sly EOS 0, H{| and polytropic EOS: 
p = lOO/? 1 - 8 , respectively. It is clearly seen that (|7.2I) indeed leads to accurate numerical results for the p-modes 
irrespective of the EOS and the compactness of the star. Again the accuracy improves with increasing frequency, 
reminding us of the insignificance of the coupling between spacetime and matter in high order p-modes. As for the 
/-modes, we note that the performance of CA (|7.2p improves for polytropic stars with low compactness. 

For the purpose of comparison, approximate eigenfrequencies a n obtained from (|7.1I) are also listed in Tables HI ITT1 
and Mil While Eq. (|7.1I) in general can yield almost exact numerical results for p-modes because it has properly taken 
the influence of F on H into consideration, it completely fails to locate approximate position of the /-mode. The 
inability of (I7.1[) to handle the /-mode can be understood as a direct consequence of the fact that D could vanish at 
certain positions inside the star at frequencies close to the /-mode frequency (see Fig. [5]). 



VIII. INVERSE-COWLING APPROXIMATION FOR w-MODE 



A. Spacetime oscillations 

From the nature of polar u>-mode oscillations, we know that the coupling between fluid motion and spacetime 
oscillations is very weak (see, e.g., [32[). Therefore, fluid motion is hardly involved in w-mode and accordingly the 
contribution of the H and dH / dr* terms in (|6.1[) is negligible. Hence, the polar w-mode system can be described by 
a single second-order equation: 

~d^ ~ V "° aSl dr^ + (Uj2 K ° " VfS ° aS ° )S = °' (8 - 1} 

where V s q, Vj , a s \ and a s o take the forms (|6.3[) . (|6.4p . (|5.14p and (|5.15l) inside the star, and take the forms (|6.2ip . 
(16.221). (|5.24l) and (|5.25p outside the star. As the fluid variable H is eliminated, the method is another form of ICA 
[3l[ |32| . To our knowledge, this is the first time that a single second-order equation for polar w-modes has been 
proposed and verified numerically (see the following discussion). 

The eigenvalue system can be simplified a lot after turning into the S representation. The spacetime oscillations 
are described by the following equation in S representation, 

S(n)=0, (8.2) 
where V takes the form V' 1 ' (16. 14)) and T^( c ) (|6.24p inside and outside the star, respectively. 



— +w 2 - V(w,r„) 



B. Numerical results and discussions 



The spacetime w-modes can be located by solving (|8.ip or equivalently (|8.2I) with the outgoing boundary condition. 
The numerical results of ICA proposed here (the columns labelled with ICA) are compared with the exact polar 
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w-mode eigenfrequencies in Table IIV1 where the exact and approximate QNM frequencies of the leading twelve polar 
w-modes (including two wn modes) for the above-mentioned EOS A NS (C = 0.27) are listed. The agreement between 
the exact frequencies and the ICA frequencies is almost perfect, especially for w r and higher order modes. 

C. High-frequency ICA 

We note that both a s o and a s i go as 1/lu 2 at high frequencies. Thus, we expect that w-modc QNMs with high 
frequencies can be approximately located by omitting terms proportional to a s o and a s ± in (|8.1I) . namely, 

uj 2 S + -| + ^ [27rr 3 (p + 3p) +m-(n + l)r] S = 0. (8.3) 

CLT ^ T 

Such a scheme is termed as the high-frequency inverse- Cowling approximation (HFICA) in the present paper. As 
shown in Table IIV1 HFICA works well in the high frequency regime. In fact, HFICA can successfully explain the 
high-frequency asymptotic behavior of polar w-modes of compact stars and an iterative method has be formulated to 
improve the accuracy of HFICA for the low-order modes [46( . 

IX. CONCLUSION AND DISCUSSION 

In this paper, based on the AAKS formalism for pulsations of compact stars (23j, we establish an equivalent 
description consisting of two coupled second-order equations in two variables, namely H and S, measuring relativistic 
dynamics of matter and spacetime, respectively. In addition to simplification of the theory proposed in [23| , the H-S 
approach developed here readily leads to CA (ICA) with unprecedented accuracies for polar p-modes (w-modes) when 
the two second-order equations are decoupled. Besides, under the CA the wave functions H n (n = 0,1,2,.. .) of the 
fluid modes (including the /-mode and the p-modes) are shown to form a complete orthonormal set, which paves the 
way for further developing perturbation schemes to improve the accuracies of the CA and to include the effect of GW 
radiation damping. 

We note that there are other theories for pulsations of compact stars which also result in two coupled second- 
order equations. For example, in the LMI formalism [2l[ the metric variable Ho(r,t) and the fluid variable 8U(r,t) = 
5p(r, t)/(p+p) + Ho(r, t)/2 are used to formulate a relativistic theory for stellar pulsations. Similar to the development 
of the present paper, the two coupled second-order equations in the LMI formalism can be decoupled by adopting the 
approximations Ho(r, t) = and 8U = 0, leading to another version of CA and ICA, respectively. The CA developed 
in the LMI formalism was applied to find the frequency of the /-mode of a polytropic NS with M = 1.4M Q and the 
percentage error was about 20% for the case I = 2 [21j . which is similar to ours. We also compute the QNM frequencies 
of the leading fluid modes of the EOS A star considered above (C = 0.27) with the CA of the LMI formalism and 
list them in Table |U It is seen that the CA in the LMI formalism can yield accurate frequencies for the p-modes. On 
the other hand, the frequencies of the lea ding polar w-modes obtained from the ICA scheme suggested in the LMI 
formalism (the columns labelled with Ref. [2l|) are listed in the Table llVl It is clear, from the numerical results, that 
the performance of the ICA scheme mentioned in [5l[ is far from satisfactory. 

The success of the CA and ICA established in the present paper based on the AAKS formalism [23| is attributable 
to the selection of independent variables, namely H and S. While the former is proportional to the change in the 
Eulerian pressure, i.e., a proper Newtonian quantity, the latter vanishes in the weak field limit. Thus, the interplay 
between H and S in fluid modes and spacetime modes is expected to be negligible, which is numerically verified in the 
above discussion. According to our experience, the choice of the independent variables H and S is the most crucial 
factor lea ding to the high precision of the ICA scheme developed here. We note that several other schemes for CA 
(see, e.g., [30(| and references therein) and ICA [3l[ based on the LD formalism have been proposed. Yet, the accuracy 
of these schemes is worse than that achieved in the present paper. 

Lastly, as F is proportional to the change in the Newtonian potential in the weak field limit, it can also reveal the 
motion of matter indirectly and replace the role of H . In fact, one can also formulate another equivalent description 
for stellar pulsation using the variables F and S to develop the F-S approach based on the AAKS formalism (see 
Appendix [51 for detailed discussion). Similar to the case in the LMI formalism mentioned above [2l|, such an F-S 
approach can yield valid CA with accuracies worse than those obtained from (|7.2[) and (|7.ip . but fails to predict 
accurate polar w-modes. The ICA in the F-S approach is identical to the HFICA in the H-S approach and hence 
works well only in the high frequency limit (see Appendix [B] and Table IIV|) . Instead, according to our experience, 
Eq. (I8.ip . or equivalently (I8.2[) . can locate all polar w-modes including the low- frequency wn modes with excellent 
numerical accuracies. In a nutshell, the H-S approach proposed here based on the AAKS formalism [23[ successfully 
leads to feasible and accurate CA and ICA schemes in a unified framework. 
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Appendix A: Boundary conditions at the stellar surface 



Near the stellar surface, we can assume that the behavior of the matter component is well captured by the polytropic 
EOS p — Kp 1 , where n > is a constant and 7 = 1 + 1/N, with N being the effective polytropic index there. Without 
loss of generality and also due to the behavior of nuclear matter at low densities, we assume that N > 1. It follows 
directly from the TOV equations that [I?} 

(Al) 
(A2) 



/-N+l 1 /-N+2 

p = aQ + a\Q 



P = /?C Ar + /?iC iV+i + 



JV+1 



with £ = R — r and 



M 



N+l 



Oil 

Pi 



kR(R-2M)(N +1) 
[M - 2{N + 1)(M - R)]M N+1 
2k n [R{R - 2M)] N + 2 {N + 1) N + 1 

M " ^ 



kR(R~ 2M)(N + 1) 
N[M - 2(N + 1)(M - R)}M N 
2k n [R{R-2M){N + 1)} n+1 ' 



Similarly, it can also be shown that 

m(r) = 



f ^r' 2 p(r')dr' = M - 4n(3R 2 C 
Jo 



2aJV+1 



2M 
~R~ 



2M( 



R 



1 



R - 2M 
M 



R(R - 2M) 
2M( 



R(R-2M)_ 

M[M -2{N + \){M - R)] 



-C + 

R(R-2M)N* 2N(N + l)R 2 (R- 2Mf s 



where 



NR 2 

C* 



37V 



-x + 

NR 



N 



2(N + 1) 



exp{[A(r') - v(r')] /2}dr' 



R-2M 



1 



2M( 



(A3) 
(A4) 
(A5) 
(A6) 

(A7) 
(A8) 
(A9) 
(A10) 
(All) 



(A12) 



R(R-2M) 

Now we consider the behavior of the wave functions H, S and F near the stellar surface. Since the coefficients in 
evolution equations (e.g., p, p and m) are expanded in powers of C and ( N , the wave functions near the surface are 
accordingly given by: 



m,n— 

cr _ i-bfb'N fmynN 
m,n— 

rp fCfc'N \ ^ r >-m, j-nN 

m,n— 



(A13) 
(A14) 
(A15) 
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where a, a', b, c, d, h mn , s mn and f mn are constants to be determined from the AAKS equations and the 
boundary conditions at the stellar surface, and h Q0 , soo and /oo are non-zero by definition. Based on the boundedness 
of H(r — R) and the AAKS equations, we can show that all the powers a, a', b, b' , c, and d vanish. 

In principle all the coefficients h mn , s mn and f mn can be determined from the modified Hamiltonian constraint 
(|4.14|) and the AAKS equations (|4.T[) and (14. 8 [I . However, for the sake of the present paper, we only need to find the 
relevant expansions to leading few orders as shown below: 

S w -soo + sioC* + S20C* + soiCf + snC +1 , (A16) 
^ « /oo + AoC* + /20C* + /oiCf + /iiCf +1 , (A17) 
H Pa fcoo + frioC*, (A18) 



with 



i?V + AMRW - M 2 R(2 + 2n + 7i? 2 w 2 ) + 4M 3 (ti + 1) 

510 " RS[RW-M(n+1)] /0 ° (A19) 

(n + 1)R 4 uj 2 + MR(R 2 u 2 - 1 - n) + 2M 2 (n + 1) 

i? 2 [i? 3 W 2 -M(n + l)] s oo, (A20) 

soi = 0, (A21) 

s n = 0, (A22) 

If 2 2(-R-2M) r „, , 2M(E-2M)(7M-3E) , , A „„, 
S20 = — |w +^-^4 '-[M-(n + l)R]js o > L foo, (A23) 

M{n+l)(2M - R) + R 3 [M - {n + 2)R]uj 2 
■' 10 ~ R 2 [R 3 uj 2 - M(n + 1)] 

(n+l)i?(2M + ni?) + i^ 2 
fc^ur — M (n + 1) 

/01 - 0, (A25) 



{R-2M) N+1 
R N - 2 M(N + 1)' 



fn = -W ± N _ 2 > , „ ftoo, (A26) 



/20 = 4 I- 2 + [M - (n + m\ /oo + ^^00, (A27) 



2\ i? 4 L y J J JUU R 

R R 2 2M + R , 1 , #^ 2 , 

^ " 2(i? - 2M) 500 " 2(i? - 2M) 510 + ^i^ /o ° + 2i? /l0 ~ ""m"^ 00 



?7, ,4 



2M(n + 1)(2M - i?)(M + i?) + i? 3 [-13M 2 + (2n + 9)A/i? - (n + l)i? 2 ]cj 2 + i? 7 

2{R - 2M)R 3 [M(n + 1) - i? 3 w 2 ] 700 
(n + l)[6M 2 + 2(n-2)A/ J R-ni? 2 ]+i? 3 [3M + (n + l) J R]a; 2 i? 2 w 2 
+ 2(i?-2M)[M(n+l)- J R 3 W 2 ] S °° ~ l^ 00 ' (A28) 

While the expressions for soi (/01), sn (/11), and S20 (/20) are obtained by comparing the coefficients of Q^ 2 , C^ -1 
and terms, respectively, in Eq. (|4.7j) (Eq. (|4.8j) ). those for sio, /10, /iio can be found from the coefficients of 
and C^" 1 " 1 terms in the modified Hamiltonian constraint (|4.14j) . On the other hand, the values of soo, /001 ^00 can be 
considered as the three degrees of freedom of the AAKS equations and other expansion coefficients are all expressible 
in terms of them. Together with the two degrees of freedom arising at the origin, there are five degrees of freedom for 
the solution of the AAKS equations. 

With the expansions (|A16I) and (IA17[) obtained above, we can match them with the ones obtained from numerically 
integrating (14.71) and (|4.8I) outward from the origin at a point near the surface of the star where the expansions are 
valid. The continuity of S, dS/dr*, F and dF/dr* lead to four independent equations, thus reducing the degree of 
freedom of the solution from five to one. As a result, the solutions for S, F and H (evaluated from the modified 
Hamiltonian constraint) are determined up to a common multiplicative constant. 
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On the other hand, there is another group of solution with divergent fluid variable H at the stellar surface: 



m,n— 

r / . ,/mnS* S># ; 

m,n— 

H = C N+1 E hmnCTC 



m,n—0 

To the leading few orders, the wave functions are given by: 

S = SQO + hoC* + S20C*, 

F = foo + AoC* + /20C* i 
H = ZiooC* N+1 , 

with 

sio = sio, 
soi = soi = 0, 
S20 = s 2 o, 
sn = sn = 0, 

flO = JW) 

foi = foi = 0, 



2(1 



, AR-2M) N+1 ~ 
h ° ~ * R N ~ 2 M °°- 



(A29) 
(A30) 
(A31) 



(A32) 
(A33) 
(A34) 

(A35) 
(A36) 
(A37) 
(A38) 
(A39) 
(A40) 

(A41) 



Appendix B: F-S Approach 
1. Formalism 



As both F and H could reveal the Newtonian aspect of stellar pulsations, it is natural to expect that an alternative 
approach using F and S to study stellar pulsations is possible. In fact, it is straightforward to make use of the 
time-independent Hamiltonian constraint (|4.13|) to eliminate the matter field H(r*) in (|4.8|) [39j: 



w 2 F d 2 F f 1 

~cT + ~d^[ ~ \ ~ c 2 



dF 



— (m + iirpr , 

r z dr* 



"' f 1 ~(p ) -2(»+l)r 



1 



1 



,-(f+A)/2 



as 



2e 



-x 



1 - ) (« + l)-8vr(p + p)r 2 



F 



(Bl) 



(|4.7j) and (|Blj) together form two coupled second-order equations in F and S, which can be used to locate QNMs of 
compact stars and lead to the so-called F-S approach in the present paper. 

Similar to the H-S approach, the S terms in (|B1[) are expected to be negligible for fluid mode pulsations, yielding 



■F d 2 F 



1 



,lr- V C 2 



f>+*)/ 2 , 3N dF 

- (to + Anpr' ) - — 

r z dr* 



I •"-/' ' pr) m\l --21" + I)/' 



F = 0. 



(B2) 



When subjected to suitable boundary condition (as discussed in the following), (|B2|) could give good approximation 
to p-mode pulsations of realistic compact stars, thereby resulting in the CA for p-modes in the F-S approach. 

On the other hand, we could also obtain the ICA for polar w-modes within the F-S approach by neglecting the F 
term in (|4.7p . i.e., 



d 2 S 2e v 



Ul 



S + -i^ + — [27rr 3 (p + 3p) + m - (n + l)r] 5" = 0. 



(B3) 
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Interestingly, such an ICA is merely the HFICA in the H-S approach, which works well only for w-modes with high 
frequencies (see Table HVl for the accuracies and validity of HFICA in the H-S approach). It reflects the fact that the 
metric variable F also participates in the dynamics of stellar pulsations and its contribution to (|4.7|) is non-negligible 
for w-modes with low frequencies. 

2. Cowling approximation 

Within the F-S approach, p-mode pulsations of compact stars in the CA can be located by solving (|B2|) with 
properly imposed boundary conditions at the two endpoints r — and r = R. As usual, the regularity boundary 
condition F oc r\ holds around the origin. On the other hand, by assuming the polytropic EOS mentioned in Appendix 
IA1 near the stellar surface, Eq. (|B2[) becomes 



d 2 F (N \ dF 



— + C- 



F = 0, 



with 



.4 

B 

C 



M 
R 2 

,2 E>2 



N + l 



NR 



N 



2(N+1) 



lo 2 R 2 N 3N(R - 2M) 
M + R 2 ' 
NR N 
~W ~ 2{N + 1) 



3N 



R-2M 
i? 4 



M(l - 9N) 



3NM 
2{N+1) 



2(n + l)R-6NR 



Correspondingly, there are two possible solutions near the surface, namely F a ~ fo + fiC* + /2C* 
/nCf +1 + ■ • where 

_ u 2 R 2 3CR-2M) 
h - [— rr- + m J/o, 



M 
fo 



R 2 
2M - R 



[M(3N - 1) - R(l + 7N)} - ^[M(8N - 1) - 4NR] + 



N -1 \ R 4 

In the CA, S is neglected and hence the modified Hamiltonian constraint (|4. 14|) leads to 



iirre u (p + p)H — 



e (f+A)/2 dF e u 

(m + 47rr 3 p)- \-uj 2 F+ — Unr 3 p + 3m) F. 



(B4) 

(B5) 

(B6) 
(B7) 
and Fb ~ 

(B8) 
(B9) 

(BIO) 



By inserting the asymptotic forms F a and Fb into the above equation, the corresponding expansion of H at the stellar 
surface can be obtained, namely, H a ~ /10C* ^ N+1 ^ and Hb ~ hoo, where 

2(2M - R) 



i? 4 



[2NM - (n + 2)iZ + 7V(n - 6)i?] 



-(4i?-7M) 



ilf 



M 2 



fo 



Q (R-2M) N + 1 N(N-1) , 
" 87r P TTW^l ~ h °> 



R 



M 



f 
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(R-2M)^ 



(Bll) 
(B12) 



As the fluid variable H is bounded, the correct boundary condition near the stellar surface is F ~ /nCi V+1 + ' ' '• 

Upon imposing the two boundary conditions mentioned above on the solution of (|F32|) , we can find the approximate 
values of the real part of the eigenfrequencies of p-modes. However, as will be discussed and explained in the following, 
we cannot locate the /-mode oscillation under the CA in this approach. 



3. Completeness and orthogonality of p-modes 



Starting from CA (IB2j) in the F-S approach, we can get a standard Sturm-Liouville equation in terms of the variable 

G = F/p, 



d 



Pg (r*)-T- 

ctr* 



-Q G (n)G + w 2 A G (r,)G = 0, 



(B13) 
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with 



Qo(r*) 
A G (r*) 



e»(p + pY 



2r(r — 2m) 



[l67rmr(p + 3p) - 647r 2 r 4 p 2 - 8irr 2 (p + p) + l] + 



(r — 2m)p 



P 



r 3 (p + p) 

„2 



47rr 3 (3p + -£) - m(l - ^) - 2(n + 1) 



C 2 



C 2 ' 



p 



e»(p + p)Cf 



(B14) 

(B15) 
(B16) 



We note that the coefficients Pair*) and Kg{t*) are strictly positive, as required for a normal Sturm-Liouville 
equation (see, e.g., [Hj]). Besides, Aq vanishes at both r = and r = R. Following directly from the physical 
boundary condition on F, namely, F ~ C^" 1 " 1 near the stellar surface, G tends to a finite constant there. In fact, by 
expanding (|B13|) about the regular singular point at r* = i?*, we find the corresponding boundary condition for G: 



l_dG_ 

G dr* 



2(4n+l)Af-(5n + 2)i? nR 2 u? 



r,=R, 



(n + 2)R 2 



(n + 2)M' 



(B17) 



Thus, a standard Sturm-Liouville eigenvalue system is established and the normalized eigenfunctions of (|B13p . G n (r*), 
where n = 1,2,3,..., form a complete set and satisfy the orthogonality relation (see, e.g., }43j): 



p 2 G m {r*)G n {r*) 
°S(p + p)C 2 



(B18) 



4. Numerical results 



The approximate real eigenfrequencies a n obtained from the CA scheme (|B2I) are listed and compared with the exact 
values uj t in Tables III \H\ and Hill For p-modes of NSs constructed with realistic EOSs, e.g., A and SLy (see Tables U 
and |lTl , the percentage error in tu r is usually less than a few percent and decreases with the mode order and the 
compactness of the star. However, the accuracy of the CA scheme (jB2[) could substantially worsen for polytropic stars 
(especially those with large polytropic indices) in spite of the fact that it still improves with increasing compactness 
and mode order (see Table ITO]) . In comparison with the two CA schemes proposed in the H-S approach, i.e., (|7.2j) 
and (|7.1[) . (jB2[) is in general the least accurate one. 

Moreover, it is worthwhile to note that the CA scheme (jB2[) fails to locate the /-mode. Hence, the n-th eigenfunction 
G n (n = 1,2,3,...) in fact corresponds to the n-th p-mode. In Fig. [3] we show the scaled eigenfunction G n defined as: 

G n ( u ) = pGn{r * ] , (B19) 
^e"(p+p)Cr 

which satisfies the orthonormal condition J^* G m (r*)G n (r*)drt, — <5 mn . There are n — 1 nodes (excluding the two 
endpoints r — and r = R) in the wave function G n or equivalently G n . In particular, the wave function of the 
pi-mode is nodeless. This observation further confirms that the pi-mode is indeed the ground state of the present 
eigen-system (see, e.g., (43|) and explains why the /-mode is missing in the CA scheme (jB2l) of the F-S approach. 
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TABLE I: Exact QNM frequency (w r , cji) and the CA frequencies a„ 
for the leading fluid modes of a EOS A star (C = 0.27). 



obtained respectively from (|7.2|) and Ref. [21J are compared 



Mode 



Exact 
(w r , Wi) 



[73] [7X) (B2) Rcf. [21] 



/ (0.1447, 7.541E-5) 0.1163 NA NA 0.1665 

Pl (0.3932, 3.047E-6) 0.3928 0.3934 0.3903 0.3984 

p 2 (0.6130, 5.824E-7) 0.6135 0.6130 0.6126 0.6159 

p 3 (0.7639, 5.269E-8) 0.7640 0.7639 0.7636 0.7664 



TABLE II: Exact complex eigenfrequencies (ui r ,uji) and approximate real eigenfrequencies a„ obtained respectively from CA 

(EU and (|B2| are listed for the leading fluid modes of NSs (EOS: Sly [HE!) 
labelled by 'NA' signify modes that cannot be located by a specific approximation. 



with different compactness C. Entries 



C Mode 



0.117 / 

0.117 pi 

0.117 p 2 

0.117 p 3 



0.210 / 

0.210 pi 

0.210 p 2 

0.210 p 3 



0.276 / 

0.276 pi 

0.276 p 2 

0.276 p 3 



Exact 

(cJr,CJi) 



l7~2l ED CHUT 



(0.0488 
(0.1654 
(0.2079 
(0.2390 



(0.1047 
(0.3281 
(0.4995 
(0.6154 



(0.1476 
(0.4074 
(0.6236 
(0.8192 



1.162E-5 
4.866E-7 
2.029E-8 
5.709E-9 



5.290E-5 
2.371E-6 
3.435E-8 
6.028E-9 



7.570E-5 
1.884E-6 
6.806E-7 
7.455E-8 



0.0410 
0.1651 
0.2079 
0.2390 



NA 
0.1654 
0.2079 
0.2390 



NA 
0.1628 
0.2071 
0.2383 



0.0848 
0.3279 
0.4996 
0.6155 



NA 
0.3283 
0.4995 
0.6154 



NA 
0.3245 
0.4985 
0.6150 



0.1177 
0.4072 
0.6241 
0.8195 



NA 
0.4075 
0.6236 
0.8192 



NA 
0.4053 
0.6234 
0.8193 



TABLE III: Exact complex eigenfrequencies (u r ,cji) and approximate real eigenfrequencies a„ obtained respectively from CA 

with different compactness C. 



(|7.2p . (|7.1|) and (|B2|) are listed for the leading fluid modes of polytropic stars (EOS: p = lOOp 
Entries labelled by 'NA' signify modes that cannot be located by a specific approximation. 



C Mode 



0.117 / 

0.117 pi 

0.117 p 2 

0.117 p 3 



0.204 / 

0.204 pi 

0.204 p 2 

0.204 p 3 



0.271 / 

0.271 pi 

0.271 p 2 

0.271 p 3 



Exact 
(w r , U)i) 



HU) (ED p2f 



(0.0525, 
(0.1160, 
(0.1761, 
(0.2338, 



1.329E-5 
3.283E-6 
2.680E-7 
1.473E-8 



(0.1160, 
(0.2337, 
(0.3477, 
(0.4584, 



5.771E-5 
1.728E-5 
5.490E-7 
3.069E-8 



(0.1659, 
(0.3012, 
(0.4249, 
(0.5496, 



5.795E-5 
1.341E-5 
2.832E-5 
1.359E-5 



0.0475 
0.1149 
0.1757 
0.2336 



NA 
0.1164 
0.1762 
0.2338 



NA 
0.1074 
0.1727 
0.2320 



0.1021 
0.2307 
0.3470 
0.4582 



NA 
0.2343 
0.3478 
0.4584 



NA 
0.2202 
0.3429 
0.4561 



0.1473 
0.2948 
0.4242 
0.5498 



NA 
0.3020 
0.4250 
0.5495 



NA 
0.2925 
0.4249 
0.5501 
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TABLE IV: QNM frequencies obtained from various ICA schemes for the leading polar ui-modes of a EOS A star (C = 0.27). 
The first two modes are classified as the wn modes. Entries labelled by 'NA' signify modes that cannot be located by a specific 
approximation. 

Exact Exact ICA ICA HFICA HFICA Rcf. [21] Ref. [21] 



CJr UJi UJi UJi UJi UJi UJi UJi 






.2067 


0. 


.7935 





.2065 





.7932 





.3304 





.8922 


NA 


NA 





.4593 





.3891 





.4599 





.3884 





.5076 





.4980 


NA 


NA 





.5096 





.1777 





.5095 





.1781 





.5102 





.1131 


0.4348 


0.0605 





.8730 





.3025 


0. 


.8730 





.3028 


0. 


.8335 





.2781 


0.7924 


0.1244 


1. 


.2266 





.3507 


1. 


.2266 





.3508 


1. 


.2006 





.3416 


1.1446 


0.1555 


1. 


.5782 





.3837 


1. 


.5782 





.3838 


1. 


.5589 





.3788 


1.4961 


0.1749 


1 


.9302 





.4093 


1 


.9301 





.4093 


1 


.9148 





.4062 


1.8476 


0.1892 


2 


.2828 





.4308 


2 


.2828 





.4308 


2 


.2700 





.4287 


2.2004 


0.2013 


2. 


.6359 





.4500 


2 


.6359 





.4500 


2 


.6250 





.4485 


2.5537 


0.2122 


2. 


.9894 





.4682 


2 


.9894 





.4683 


2 


.9799 





.4670 


2.9079 


0.2231 


3 


.3429 





.4860 


3 


.3428 





.4860 


3 


.3344 





.4850 


3.2619 


0.2347 


3 


.6960 





.5038 


3 


.6960 





.5038 


3 


.6884 





.5030 


3.6158 


0.2468 
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FIG. 1: The normalized //-functions for (a) the /-mode; (b) the pi-mode; and (c) the p2-mode of a star constructed with 
EOS A (p c = 2.227 x 10 -3 km -2 , C = 0.27) are plotted against the normalized radius r/R. The solid and the dashed lines are 
respectively the real and the imaginary parts of //-function under CA (|7.2[1 . The dots and crosses respectively show the real 
and the imaginary parts of the exact //-function. 
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FIG. 2: The function D(r,uj) is plotted against r/R with uj — <7o,cri,<72. 




FIG. 3: The scaled function G„ of a NS with SLy EOS 0, 5§] and compactness C = 0.276 is plotted against r/R for n = 1 
(solid line), 2 (dashed line), 3 (dot-dashed line), which correspond to pi, P2 and p3 modes, respectively. 



